Exploring the mechanisms underlying stroke volume variability reduction in a murine model of heart failure with reduced ejection fraction.

Heart failure with reduced ejection fraction (HFrEF) is accompanied by disregulation of cardiovascular function. Heart rate variability (HRV) is commonly used to assess autonomic dysfunction in HFrEF. However, analysis of stroke volume variability (SVV) may provide additional insights. We examined HRV and SVV in a mouse model of HFrEF. HFrEF mice exhibited reduced stroke volume and ejection fraction versus controls, confirming cardiac contractile dysfunction. HRV was preserved in HFrEF mice. However, SVV was markedly diminished, indicating dissociation between HRV and SVV regulation. Using a mathematical model, we propose that Frank-Starling mechanism abnormalities in HFrEF disrupt SVV independent of HRV. Assessing SVV could thus provide unique insights beyond HRV into cardiovascular control deficits in HFrEF.

Thermodynamic performance of coupled enzymatic reactions: A chemical kinetics model for analyzing cotransporters, ion pumps, and ATP synthases.

Previous research has suggested that molecular energy converters such as ATP synthases, ion pumps, and cotransporters operate via spatially separate pathways for free energy donor and acceptor reactions linked by a protein molecule. We present a chemical kinetics model based on these works, with the basic assumption that all molecular energy converters can be thought of as linked enzymatic reactions, one running downhill the chemical potential gradient and driving the other uphill. To develop the model we first look at how an enzyme process can be forced to go backwards using a basic kinetic model. We then use these findings to suggest a thermodynamically consistent method of linking two enzymatic reactions. Finally, in the context of the aforementioned energy converters, the thermodynamic performance of the resulting model is thoroughly investigated and the obtained results are contrasted with experimental data.

On the effects of diabetes mellitus on the mechanical properties of DRG sensory neurons and their possible relation with diabetic neuropathy.

Diabetic neuropathy (DN) is one of the principal complications of diabetes mellitus (DM). Dorsal root ganglion (DRG) neurons are the primary sensory neurons that transduce mechanical, chemical, thermal, and pain stimuli. Diabetes-caused sensitivity alterations and presence of pain are due to cellular damage originated by persistent hyperglycemia, microvascular insufficiency, and oxidative and nitrosative stress. However, the underlying mechanisms have not been fully clarified. The present work addresses this problem by hypothesizing that sensitivity changes in DN result from mechanotransduction-system alterations in sensory neurons; especially, plasma membrane affectations. This hypothesis is tackled by means of elastic-deformation experiments performed on DGR neurons from a murine model for type-1 DM, as well a mathematical model of the cell mechanical structure. The obtained results suggest that the plasma-membrane fluidity of DRG sensory neurons is modified by the induction of DM, and that this alteration may correlate with changes in the cell calcium transient that results from mechanical stimuli.

Hybrid reaction-diffusion and clock-and-wavefront model for the arrest of oscillations in the somitogenesis segmentation clock.

The clock and wavefront paradigm is arguably the most widely accepted model for explaining the embryonic process of somitogenesis. According to this model, somitogenesis is based upon the interaction between a genetic oscillator, known as segmentation clock, and a differentiation wavefront, which provides the positional information indicating where each pair of somites is formed. Shortly after the clock and wavefront paradigm was introduced, Meinhardt presented a conceptually different mathematical model for morphogenesis in general, and somitogenesis in particular. Recently, Cotterell et al. [A local, self-organizing reaction-diffusion model can explain somite patterning in embryos, Cell Syst. 1, 257-269 (2015)] rediscovered an equivalent model by systematically enumerating and studying small networks performing segmentation. Cotterell et al. called it a progressive oscillatory reaction-diffusion (PORD) model. In the Meinhardt-PORD model, somitogenesis is driven by short-range interactions and the posterior movement of the front is a local, emergent phenomenon, which is not controlled by global positional information. With this model, it is possible to explain some experimental observations that are incompatible with the clock and wavefront model. However, the Meinhardt-PORD model has some important disadvantages of its own. Namely, it is quite sensitive to fluctuations and depends on very specific initial conditions (which are not biologically realistic). In this work, we propose an equivalent Meinhardt-PORD model and then amend it to couple it with a wavefront consisting of a receding morphogen gradient. By doing so, we get a hybrid model between the Meinhardt-PORD and the clock-and-wavefront ones, which overcomes most of the deficiencies of the two originating models.

Dynamics of Mammalian Cell Infection by Trypanosoma cruzi trypomastigotes.

In a recent work we demonstrated that Trypanosoma cruzi trypomastigotes change their motility patterns in the presence of mammalian cells, that the extent of the changes depends on the cell line, and that this extent is positively correlated with the efficiency with which parasites invade the different cell lines. These results open the question of what cellular characteristics are relevant for parasite identification and invasion. In the present work, we tackled such question. We performed infection-kinetics experiments on various cell lines, and developed a mathematical model to simulate the experimental outcomes. An analysis of the cell-parasite mechanisms included in the model, together with the parameter values that allowed it to replicate the experimental results, suggests that a process related to the cell replication rate may strongly influence the parasite invasion efficiency, and the infection dynamics in general.

Motility patterns of Trypanosoma cruzi trypomastigotes correlate with the efficiency of parasite invasion in vitro.

Numerous works have demonstrated that trypanosomatid motility is relevant for parasite replication and sensitivity. Nonetheless, although some findings indirectly suggest that motility also plays an important role during infection, this has not been extensively investigated. This work is aimed at partially filling this void for the case of Trypanosoma cruzi. After recording swimming T. cruzi trypomastigotes (CL Brener strain) and recovering their individual trajectories, we statistically analyzed parasite motility patterns. We did this with parasites that swim alone or above monolayer cultures of different cell lines. Our results indicate that T. cruzi trypomastigotes change their motility patterns when they are in the presence of mammalian cells, in a cell-line dependent manner. We further performed infection experiments in which each of the mammalian cell cultures were incubated for 2 h together with trypomastigotes, and measured the corresponding invasion efficiency. Not only this parameter varied from cell line to cell line, but it resulted to be positively correlated with the corresponding intensity of the motility pattern changes. Together, these results suggest that T. cruzi trypomastigotes are capable of sensing the presence of mammalian cells and of changing their motility patterns accordingly, and that this might increase their invasion efficiency.

Crosstalk dynamics between the circadian clock and the mTORC1 pathway.

Crosstalk between the circadian clock clockwork and cellular metabolic regulatory networks is crucial to ensure an adequate response of an organism to the day/night cycle. mTOR (mammalian/mechanistic target of rapamycin) is a master growth regulator and sensor of nutrient status, which is part of the mTOR complex 1 (mTORC1). While the circadian clock confers rhythmicity to the mTOR protein by regulating its degradation rate, mTORC1 activity diminishes period and augments amplitude of circadian oscillations at the cellular level by a currently unknown mechanism. Here, we develop a mathematical deterministic DAE (differential-algebraic equation) model, to explore the possible interactions that allow mTORC1 to display such regulation of the core circadian clock. Our results suggest that mTORC1 is capable of regulating amplitude by exerting translational control on core the clock protein BMAL1, and that period-tuning is achieved by controlling post-translational localization of BMAL1. Since, in our model, mTORC1 control of BMAL1 localization greatly diminishes the ability of the clock to oscillate, and regulation of BMAL1 translation reduces this effect, our results also suggest that both levels of regulation must be present to ensure the robustness of oscillations. Together, the above results emphasize the importance of the influence of mTORC1 on the circadian rhythms.

Three-Way Interactions in an Artificial Community of Bacterial Strains Directly Isolated From the Environment and Their Effect on the System Population Dynamics.

This work is motivated by previous studies that have analyzed the population ecology of a collection of culturable thermoresistant bacteria, isolated from the Churince lagoon in Cuatro Cienegas, Mexico. In particular, it is aimed at testing a hypothesis from a modeling study, which states that antagonistic and sensitive bacteria co-exist thanks to resistant bacteria that protect sensitive ones by forming physical barriers. We selected three different bacterial strains from the referred collection: one antagonistic, one sensitive, and one resistant, and studied the population dynamics of mixed colonies. Our results show that, although the proposed protective mechanism does not work in this case, the resistant strain confers some kind of protection to sensitive bacteria. Further modeling and experimental results suggest that the presence of resistant bacteria indirectly improves the probability that patches of sensitive bacteria grow in a mixed colony. More precisely, our results suggest that by making antagonistic bacteria produce and secrete an antagonistic substance (with the concomitant metabolic cost and growth rate reduction), resistant bacteria increase the likelihood that sensitive bacteria locally outcompete antagonistic ones.

Polymerization of sarcoplasmic-reticulum calcium-binding proteins might explain observed reticulum kinetics-on-demand behavior.

Reported experimental results, in which transient elevations of sarcoplasmic calcium levels are induced by caffeine in smooth muscle cells, apparently contradict the principle of mass conservation. The commonly accepted model assumes that the total number of Ca2+ binding sites is fixed. A former work dealing with this problem proved that assuming the presence within the reticulum of calcium sequestering proteins whose total number of calcium binding sites increases as the existent sites get occupied, is enough to explain the above referred counter-intuitive experimental results. However, no chemical explanation was given to account for this binding-site count increase. In the present work, we propose a chemical-kinetics scheme for the binding of calcium to calsequestrin (a protein found within the reticulum) and the polymerization of this protein. On the one hand, this scheme is in agreement with reported results on calsequestrin binding kinetics, but it is also fully capable of explaining the observed intriguing performance of the sarcoplasmic reticulum. We further explore the behavior of the resulting nonlinear dynamic system and discuss possible physiological implications of the proposed scheme.

Bistable behaviour and medium-dependent post-translational regulation of the tryptophanase operon regulatory pathway in Echerichia coli.

The present work is aimed at studying the dynamic behaviour of the tryptopnanase (tna) operon, which encodes the proteins necessary to uptake and metabolise tryptophan to use it as a carbon source in the absence of glucose. To this end, we designed a micro-bioreactor capable of driving a bacterial culture to a stationary state. This allowed us to explore (at the single cell level) the tna operon steady-state dynamics under multiple culture conditions. Our experimental results suggest that the tna operon is bistable for a specific range of environmental tryptophan and glucose concentrations, and evidence that both reagents play a role on the activation of the enzyme in charge of metabolising tryptophan: tryptophanase (TnaA). Based on our experimental data and the already known regulatory mechanisms, we developed a mathematical model for the tna operon regulatory pathway. Our modelling results reinforce the claim that the tna operon is bistable, and further suggest that the activity of enzyme TnaA is regulated by the environmental levels of glucose and tryptophan via a common signalling pathway. Possible biological implications of our findings are further discussed.

Optimization and Stability of Heat Engines: The Role of Entropy Evolution.

Local stability of maximum power and maximum compromise (Omega) operation regimes dynamic evolution for a low-dissipation heat engine is analyzed. The thermodynamic behavior of trajectories to the stationary state, after perturbing the operation regime, display a trade-off between stability, entropy production, efficiency and power output. This allows considering stability and optimization as connected pieces of a single phenomenon. Trajectories inside the basin of attraction display the smallest entropy drops. Additionally, it was found that time constraints, related with irreversible and endoreversible behaviors, influence the thermodynamic evolution of relaxation trajectories. The behavior of the evolution in terms of the symmetries of the model and the applied thermal gradients was analyzed.

Motility of Escherichia coli in a quasi-two-dimensional porous medium.

Bacterial migration through confined spaces is critical for several phenomena, such as biofilm formation, bacterial transport in soils, and bacterial therapy against cancer. In the present work, E. coli (strain K12-MG1655 WT) motility was characterized by recording and analyzing individual bacterium trajectories in a simulated quasi-two-dimensional porous medium. The porous medium was simulated by enclosing, between slide and cover slip, a bacterial-culture sample mixed with uniform 2.98-µm-diameter spherical latex particles. The porosity of the medium was controlled by changing the latex particle concentration. By statistically analyzing several trajectory parameters (instantaneous velocity, turn angle, mean squared displacement, etc.), and contrasting with the results of a random-walk model developed ad hoc, we were able to quantify the effects that different obstacle concentrations have upon bacterial motility.

Experimental and Mathematical-Modeling Characterization of Trypanosoma cruzi Epimastigote Motility.

The present work is aimed at characterizing the motility of parasite T. cruzi in its epimastigote form. To that end, we recorded the trajectories of two strains of this parasite (a wild-type strain and a stable transfected strain, which contains an ectopic copy of LYT1 gene and whose motility is known to be affected). We further extracted parasite trajectories from the recorded videos, and statistically analysed the following trajectory-step features: step length, angular change of direction, longitudinal and transverse displacements with respect to the previous step, and mean square displacement. Based on the resulting observations, we developed a mathematical model to simulate parasite trajectories. The fact that the model predictions closely match most of the experimentally observed parasite-trajectory characteristics, allows us to conclude that the model is an accurate description of T. cruzi motility.

Kinetics on Demand Is a Simple Mathematical Solution that Fits Recorded Caffeine-Induced Luminal SR Ca2+ Changes in Smooth Muscle Cells.

The process of Ca2+ release from sarcoplasmic reticulum (SR) comprises 4 phases in smooth muscle cells. Phase 1 is characterized by a large increase of the intracellular Ca2+ concentration ([Ca2+]i) with a minimal reduction of the free luminal SR [Ca2+] ([Ca2+]FSR). Importantly, active SR Ca2+ ATPases (SERCA pumps) are necessary for phase 1 to occur. This situation cannot be explained by the standard kinetics that involves a fixed amount of luminal Ca2+ binding sites. A new mathematical model was developed that assumes an increasing SR Ca2+ buffering capacity in response to an increase of the luminal SR [Ca2+] that is called Kinetics-on-Demand (KonD) model. This approach can explain both phase 1 and the refractory period associated with a recovered [Ca2+]FSR. Additionally, our data suggest that active SERCA pumps are a requisite for KonD to be functional; otherwise luminal SR Ca2+ binding proteins switch to standard kinetics. The importance of KonD Ca2+ binding properties is twofold: a more efficient Ca2+ release process and that [Ca2+]FSR and Ca2+-bound to SR proteins ([Ca2+]BSR) can be regulated separately allowing for Ca2+ release to occur (provided by Ca2+-bound to luminal Ca2+ binding proteins) without an initial reduction of the [Ca2+]FSR.

Antagonistic interactions are sufficient to explain self-assemblage of bacterial communities in a homogeneous environment: a computational modeling approach.

Most of the studies in Ecology have been devoted to analyzing the effects the environment has on individuals, populations, and communities, thus neglecting the effects of biotic interactions on the system dynamics. In the present work we study the structure of bacterial communities in the oligotrophic shallow water system of Churince, Cuatro Cienegas, Mexico. Since the physicochemical conditions of this water system are homogeneous and quite stable in time, it is an excellent candidate to study how biotic factors influence the structure of bacterial communities. In a previous study, the binary antagonistic interactions of 78 bacterial strains, isolated from Churince, were experimentally determined. We employ these data to develop a computer algorithm to simulate growth experiments in a cellular grid representing the pond. Remarkably, in our model, the dynamics of all the simulated bacterial populations is determined solely by antagonistic interactions. Our results indicate that all bacterial strains (even those that are antagonized by many other bacteria) survive in the long term, and that the underlying mechanism is the formation of bacterial community patches. Patches corresponding to less antagonistic and highly susceptible strains are consistently isolated from the highly-antagonistic bacterial colonies by patches of neutral strains. These results concur with the observed features of the bacterial community structure previously reported. Finally, we study how our findings depend on factors like initial population size, differential population growth rates, homogeneous population death rates, and enhanced bacterial diffusion.

The utility of simple mathematical models in understanding gene regulatory dynamics.

In this review, we survey work that has been carried out in the attempts of biomathematicians to understand the dynamic behaviour of simple bacterial operons starting with the initial work of the 1960's. We concentrate on the simplest of situations, discussing both repressible and inducible systems and then turning to concrete examples related to the biology of the lactose and tryptophan operons. We conclude with a brief discussion of the role of both extrinsic noise and so-called intrinsic noise in the form of translational and/or transcriptional bursting.

Optimal performance of the tryptophan operon of E. coli: a stochastic, dynamical, mathematical-modeling approach.

In this work, we develop a detailed, stochastic, dynamical model for the tryptophan operon of E. coli, and estimate all of the model parameters from reported experimental data. We further employ the model to study the system performance, considering the amount of biochemical noise in the trp level, the system rise time after a nutritional shift, and the amount of repressor molecules necessary to maintain an adequate level of repression, as indicators of the system performance regime. We demonstrate that the level of cooperativity between repressor molecules bound to the first two operators in the trp promoter affects all of the above enlisted performance characteristics. Moreover, the cooperativity level found in the wild-type bacterial strain optimizes a cost-benefit function involving low biochemical noise in the tryptophan level, short rise time after a nutritional shift, and low number of regulatory molecules.

Antagonism influences assembly of a Bacillus guild in a local community and is depicted as a food-chain network.

Understanding the principles that govern community assemblages is a central goal of ecology. There is limited experimental evidence in natural settings showing that microbial assembly in communities are influenced by antagonistic interactions. We, therefore, analyzed antagonism among bacterial isolates from a taxonomically related bacterial guild obtained from five sites in sediments from a fresh water system. We hypothesized that if antagonistic interactions acted as a shaping force of the community assembly, then the frequency of resistance to antagonism among bacterial isolates originating from a given site would be higher than the resistance to conspecifics originating from a different assemblage. Antagonism assays were conducted between 78 thermoresistant isolates, of which 72 were Bacillus spp. Sensitive, resistant and antagonistic isolates co-occurred at each site, but the within-site frequency of resistance observed was higher than that observed when assessed across-sites. We found that antagonism results from bacteriocin-like substances aimed at the exclusion of conspecifics. More than 6000 interactions were scored and described by a directed network with hierarchical structure that exhibited properties that resembled a food chain, where the different Bacillus taxonomic groups occupied specific positions. For some tested interacting pairs, the unidirectional interaction could be explained by competition that inhibited growth or completely excluded one of the pair members. This is the first report on the prevalence and specificity of Bacillus interactions in a natural setting and provides evidence for the influence of bacterial antagonist interactions in the assemblage of a taxonomically related guild in local communities.

Hierarchical dynamics of a transcription factors network in E. coli.

In the present work we study the transcription-factor regulatory network that controls the synthesis of flagella in E. coli. Our objective is to address how the transcription-factor dynamics (in terms of their promoter activities and associated rates) correlate with their positions in the hierarchical organization of this regulatory network. Our results suggest that global-regulator promoters express at higher rates than those of local regulators, particularly when the bacterial populations are actively growing. Furthermore, promoter activity decreases together with the rate of cellular division. And finally, local-regulator promoters reach their maximal activity later than global-regulator promoters do. In summary, our results suggest a strong correlation between promoter activities and their hierarchical organization in this particular regulatory network.

What is the core oscillator in the speract-activated pathway of the Strongylocentrotus purpuratus sperm flagellum?

Sperm chemotaxis has an important role in fertilization. Most of our knowledge regarding this phenomenon comes from studies in organisms whose fertilization occurs externally, like sea urchins. Sea urchin spermatozoa respond to sperm-activating peptides, which diffuse from the egg jelly coat and interact with their receptor in the flagellum, triggering several physiological responses: changes in membrane potential, intracellular pH, cyclic nucleotide levels, and intracellular Ca2+ concentration ([Ca2+]). In particular, flagellar [Ca2+] has been shown to oscillate. These [Ca2+] oscillations are correlated with changes in the flagellar shape and so with the regulation of the sperm swimming paths. In this study, we demonstrate, from a mathematical modeling perspective, that the reported speract-activated signaling pathway in Strongylocentrotus purpuratus (speract being a sperm-activating peptide specific to this species) has the necessary elements to replicate the reported [Ca2+] oscillations. We further investigate which elements of this signaling pathway constitute the core oscillator.